Geographical survey of the mycobiome and microbiome of Southern California glassy-winged sharpshooters.

The glassy-winged sharpshooter, Homalodisca vitripennis Germar, is an invasive xylem-feeding leafhopper with a devastating economic impact on California agriculture through transmission of the plant pathogen, Xylella fastidiosa. While studies have focused on X. fastidiosa or known symbionts of H. vitripennis, little work has been done at the scale of the microbiome (the bacterial community) or mycobiome (the fungal community). Here, we characterize the mycobiome and the microbiome of H. vitripennis across Southern California and explore correlations with captivity and host insecticide resistance status. Using high-throughput sequencing of the ribosomal internal transcribed spacer 1 region and the 16S rRNA gene to profile the mycobiome and microbiome, respectively, we found that while the H. vitripennis mycobiome significantly varied across Southern California, the microbiome did not. We also observed a significant difference in both the mycobiome and microbiome between captive and wild H. vitripennis. Finally, we found that the mycobiome, but not the microbiome, was correlated with insecticide resistance status in wild H. vitripennis. This study serves as a foundational look at the H. vitripennis mycobiome and microbiome across Southern California. Future work should explore the putative link between microbes and insecticide resistance status and investigate whether microbial communities should be considered in H. vitripennis management practices. IMPORTANCE The glassy-winged sharpshooter is an invasive leafhopper that feeds on the xylem of plants and transmits the devastating pathogen, Xylella fastidiosa, resulting in significant economic damage to California's agricultural system. While studies have focused on this pathogen or obligate symbionts of the glassy-winged sharpshooter, there is limited knowledge of the bacterial and fungal communities that make up its microbiome and mycobiome. To address this knowledge gap, we explored the composition of the mycobiome and the microbiome of the glassy-winged sharpshooter across Southern California and identified differences associated with geography, captivity, and host insecticide resistance status. Understanding sources of variation in the microbial communities associated with the glassy-winged sharpshooter is an important consideration for developing management strategies to control this invasive insect. This study is a first step toward understanding the role microbes may play in the glassy-winged sharpshooter's resistance to insecticides.

Genomes and transcriptomes help unravel the complex life cycle of the blastoclad fungus, Coelomomyces lativittatus, an obligate parasite of mosquitoes and microcrustaceans.

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

Highly diverse and unknown viruses may enhance Antarctic endoliths' adaptability.

BACKGROUND: Rock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities. RESULTS: We performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry. CONCLUSION: This catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate. Video Abstract.

Transcriptome and population structure of glassy-winged sharpshooters (Homalodisca vitripennis) with varying insecticide resistance in southern California.

BACKGROUND: Homalodisca vitripennis Germar, the glassy-winged sharpshooter, is an invasive insect in California and a critical threat to agriculture through its transmission of the plant pathogen, Xylella fastidiosa. Quarantine, broad-spectrum insecticides, and biological control have been used for population management of H. vitripennis since its invasion and subsequent proliferation throughout California. Recently wide-spread neonicotinoid resistance has been detected in populations of H. vitripennis in the southern portions of California's Central Valley. In order to better understand potential mechanisms of H. vitripennis neonicotinoid resistance, we performed RNA sequencing on wild-caught insecticide-resistant and relatively susceptible sharpshooters to profile their transcriptome and population structure. RESULTS: We identified 81 differentially expressed genes with higher expression in resistant individuals. The significant largest differentially expressed candidate gene linked to resistance status was a cytochrome P450 gene with similarity to CYP6A9. Furthermore, we observed an over-enrichment of GO terms representing functions supportive of roles in resistance mechanisms (cytochrome P450s, M13 peptidases, and cuticle structural proteins). Finally, we saw no evidence of broad-scale population structure, perhaps due to H. vitripennis' relatively recent introduction to California or due to the relatively small geographic scale investigated here. CONCLUSIONS: In this work, we characterized the transcriptome of insecticide-resistant and susceptible H. vitripennis and identified candidate genes that may be involved in resistance mechanisms for this species. Future work should seek to build on the transcriptome profiling performed here to confirm the role of the identified genes, particularly the cytochrome P450, in resistance in H. vitripennis. We hope this work helps aid future population management strategies for this and other species with growing insecticide resistance.

Metagenome-Assembled Genomes of Bacteria Associated with Massospora cicadina Fungal Plugs from Infected Brood VIII Periodical Cicadas.

We report six metagenome-assembled genomes (MAGs) associated with Massospora cicadina strain MCPNR19 (ARSEF 14555), an obligate entomopathogenic fungus of periodical cicadas. The MAGs include representatives of Pantoea, Pseudomonas, Lactococcus, and one potential new Chryseobacterium species. Future research is needed to resolve the ecology of these MAGs and determine whether they represent symbionts or contaminants.

An Improved 1.5-Gigabase Draft Assembly of Massospora cicadina (Zoopagomycota), an Obligate Fungal Parasite of 13- and 17-Year Cicadas.

A 1.488-Gb draft genome sequence was assembled for the fungus Massospora cicadina, an obligate parasite of periodical cicadas. The M. cicadina genome has experienced massive expansion via transposable elements (TEs), which account for 92% of the genome.

A guide to preprinting for early-career researchers.

The use of preprints, research manuscripts shared publicly before completing the traditional peer-review process, is becoming a more common practice among life science researchers. Early-career researchers (ECRs) benefit from posting preprints as they are shareable, citable, and prove productivity. However, preprinting a manuscript involves a discussion among all co-authors, and ECRs are often not the decision-makers. Therefore, ECRs may find themselves in situations where they are interested in depositing a preprint but are unsure how to approach their co-authors or advisor about preprinting. Leveraging our own experiences as ECRs, and feedback from the research community, we have constructed a guide for ECRs who are considering preprinting to enable them to take ownership over the process and to raise awareness about preprinting options. We hope that this guide helps ECRs to initiate conversations about preprinting with co-authors and encourage them to preprint their future research.

Metagenome-Assembled Genomes of Bacterial Symbionts Associated with Insecticide-Resistant and -Susceptible Individuals of the Glassy-Winged Sharpshooter (Homalodisca vitripennis).

The role of microbes in insecticide resistance is an emerging question. Here, we describe six metagenome-assembled genomes (MAGs) associated with the glassy-winged sharpshooter (Homalodisca vitripennis [Germar, 1821]) (Hemiptera, Cicadellidae). MAGs representing the obligate symbionts Candidatus Sulcia muelleri and Candidatus Baumannia cicadellinicola and the facultative symbiont Wolbachia were obtained from imidacloprid-resistant and imidacloprid-susceptible sharpshooters.

The population genetics of the causative agent of snake fungal disease indicate recent introductions to the USA.

Snake fungal disease (SFD; ophidiomycosis), caused by the pathogen Ophidiomyces ophiodiicola (Oo), has been documented in wild snakes in North America and Eurasia, and is considered an emerging disease in the eastern United States of America. However, a lack of historical disease data has made it challenging to determine whether Oo is a recent arrival to the USA or whether SFD emergence is due to other factors. Here, we examined the genomes of 82 Oo strains to determine the pathogen's history in the eastern USA. Oo strains from the USA formed a clade (Clade II) distinct from European strains (Clade I), and molecular dating indicated that these clades diverged too recently (approximately 2,000 years ago) for transcontinental dispersal of Oo to have occurred via natural snake movements across Beringia. A lack of nonrecombinant intermediates between clonal lineages in Clade II indicates that Oo has actually been introduced multiple times to North America from an unsampled source population, and molecular dating indicates that several of these introductions occurred within the last few hundred years. Molecular dating also indicated that the most common Clade II clonal lineages have expanded recently in the USA, with time of most recent common ancestor mean estimates ranging from 1985 to 2007 CE. The presence of Clade II in captive snakes worldwide demonstrates a potential mechanism of introduction and highlights that additional incursions are likely unless action is taken to reduce the risk of pathogen translocation and spillover into wild snake populations.

Efficient CRISPR/Cas9-mediated genome modification of the glassy-winged sharpshooter Homalodisca vitripennis (Germar).

CRISPR/Cas9 technology enables the extension of genetic techniques into insect pests previously refractory to genetic analysis. We report the establishment of genetic analysis in the glassy-winged sharpshooter (GWSS), Homalodisca vitripennis, which is a significant leafhopper pest of agriculture in California. We use a novel and simple approach of embryo microinjection in situ on the host plant and obtain high frequency mutagenesis, in excess of 55%, of the cinnabar and white eye pigmentation loci. Through pair matings, we obtained 100% transmission of w and cn alleles to the G3 generation and also established that both genes are located on autosomes. Our analysis of wing phenotype revealed an unexpected discovery of the participation of pteridine pigments in wing and wing-vein coloration, indicating a role for these pigments beyond eye color. We used amplicon sequencing to examine the extent of off-target mutagenesis in adults arising from injected eggs, which was found to be negligible or non-existent. Our data show that GWSS can be easily developed as a genetic model system for the Hemiptera, enabling the study of traits that contribute to the success of invasive pests and vectors of plant pathogens. This will facilitate novel genetic control strategies.

Invasive Grass Dominance over Native Forbs Is Linked to Shifts in the Bacterial Rhizosphere Microbiome.

Rhizosphere microbiomes have received growing attention in recent years for their role in plant health, stress tolerance, soil nutrition, and invasion. Still, relatively little is known about how these microbial communities are altered under plant competition, and even less about whether these shifts are tied to competitive outcomes between native and invasive plants. We investigated the structure and diversity of rhizosphere bacterial and fungal microbiomes of native annual forbs and invasive annual grasses grown in a shade-house both individually and in competition using high-throughput amplicon sequencing of the bacterial 16S rRNA gene and the fungal ITS region. We assessed how differentially abundant microbial families correlate to plant biomass under competition. We find that bacterial diversity and structure differ between native forbs and invasive grasses, but fungal diversity and structure do not. Furthermore, bacterial community structures under competition are distinct from individual bacterial community structures. We also identified five bacterial families that varied in normalized abundance between treatments and that were correlated with plant biomass under competition. We speculate that invasive grass dominance over these natives may be partially due to effects on the rhizosphere community, with changes in specific bacterial families potentially benefiting invaders at the expense of natives.

Improved draft reference genome for the Glassy-winged Sharpshooter (Homalodisca vitripennis), a vector for Pierce's disease.

Homalodisca vitripennis (Hemiptera: Cicadellidae), known as the glassy-winged sharpshooter, is a xylem feeding leafhopper and an important agricultural pest as a vector of Xylella fastidiosa, which causes Pierce's disease in grapes and a variety of other scorch diseases. The current H. vitripennis reference genome from the Baylor College of Medicine's i5k pilot project is a 1.4-Gb assembly with 110,000 scaffolds, which still has significant gaps making identification of genes difficult. To improve on this effort, we used a combination of Oxford Nanopore long-read sequencing technology combined with Illumina sequencing reads to generate a better assembly and first-pass annotation of the whole genome sequence of a wild-caught Californian (Tulare County) individual of H. vitripennis. The improved reference genome assembly for H. vitripennis is 1.93-Gb in length (21,254 scaffolds, N50 = 650 Mb, BUSCO completeness = 94.3%), with 33.06% of the genome masked as repetitive. In total, 108,762 gene models were predicted including 98,296 protein-coding genes and 10,466 tRNA genes. As an additional community resource, we identified 27 orthologous candidate genes of interest for future experimental work including phenotypic marker genes like white. Furthermore, as part of the assembly process, we generated four endosymbiont metagenome-assembled genomes, including a high-quality near complete 1.7-Mb Wolbachia sp. genome (1 scaffold, CheckM completeness = 99.4%). The improved genome assembly and annotation for H. vitripennis, curated set of candidate genes, and endosymbiont MAGs will be invaluable resources for future research of H. vitripennis.

Reconstruction of Metagenome-Assembled Genomes from Aquaria.

Here, we report 11 metagenome-assembled genomes (MAGs) reconstructed from freshwater and saltwater aquaria, including representatives of Polynucleobacter, Anaerolinea, Roseobacter, Flavobacteriia, Octadecabacter, Mycobacterium, and Candidate Phyla Radiation (CPR) members. These MAGs can serve as a resource for aquatic research and elucidating the role of CPR taxa in the built environment.

Correction: Fungi, bacteria and oomycota opportunistically isolated from the seagrass, Zostera marina.

[This corrects the article DOI: 10.1371/journal.pone.0236135.].

Global Diversity and Biogeography of the Zostera marina Mycobiome.

Seagrasses are marine flowering plants that provide critical ecosystem services in coastal environments worldwide. Marine fungi are often overlooked in microbiome and seagrass studies, despite terrestrial fungi having critical functional roles as decomposers, pathogens, or endophytes in global ecosystems. Here, we characterize the distribution of fungi associated with the seagrass Zostera marina, using leaves, roots, and rhizosphere sediment from 16 locations across its full biogeographic range. Using high-throughput sequencing of the ribosomal internal transcribed spacer (ITS) region and 18S rRNA gene, we first measured fungal community composition and diversity. We then tested hypotheses of neutral community assembly theory and the degree to which deviations suggested that amplicon sequence variants (ASVs) were plant selected or dispersal limited. Finally, we identified a core mycobiome and investigated the global distribution of differentially abundant ASVs. We found that the fungal community is significantly different between sites and that the leaf mycobiome follows a weak but significant pattern of distance decay in the Pacific Ocean. Generally, there was evidence for both deterministic and stochastic factors contributing to community assembly of the mycobiome, with most taxa assembling through stochastic processes. The Z. marina core leaf and root mycobiomes were dominated by unclassified Sordariomycetes spp., unclassified Chytridiomycota lineages (including Lobulomycetaceae spp.), unclassified Capnodiales spp., and Saccharomyces sp. It is clear from the many unclassified fungal ASVs and fungal functional guilds that knowledge of marine fungi is still rudimentary. Further studies characterizing seagrass-associated fungi are needed to understand the roles of these microorganisms generally and when associated with seagrasses. IMPORTANCE Fungi have important functional roles when associated with land plants, yet very little is known about the roles of fungi associated with marine plants, like seagrasses. In this study, we report the results of a global effort to characterize the fungi associated with the seagrass Zostera marina across its full biogeographic range. Although we defined a putative global core fungal community, it is apparent from the many fungal sequences and predicted functional guilds that had no matches to existing databases that general knowledge of seagrass-associated fungi and marine fungi is lacking. This work serves as an important foundational step toward future work investigating the functional ramifications of fungi in the marine ecosystem.

Fungi, bacteria and oomycota opportunistically isolated from the seagrass, Zostera marina.

Fungi in the marine environment are often neglected as a research topic, despite that fungi having critical roles on land as decomposers, pathogens or endophytes. Here we used culture-dependent methods to survey the fungi associated with the seagrass, Zostera marina, also obtaining bacteria and oomycete isolates in the process. A total of 108 fungi, 40 bacteria and 2 oomycetes were isolated. These isolates were then taxonomically identified using a combination of molecular and phylogenetic methods. The majority of the fungal isolates were classified as belonging to the classes Eurotiomycetes, Dothideomycetes, and Sordariomycetes. Most fungal isolates were habitat generalists like Penicillium sp. and Cladosporium sp., but we also cultured a diverse set of rare taxa including possible habitat specialists like Colletotrichum sp. which may preferentially associate with Z. marina leaf tissue. Although the bulk of bacterial isolates were identified as being from known ubiquitous marine lineages, we also obtained several Actinomycetes isolates and a Phyllobacterium sp. We identified two oomycetes, another understudied group of marine microbial eukaryotes, as Halophytophthora sp. which may be opportunistic pathogens or saprophytes of Z. marina. Overall, this study generates a culture collection of fungi which adds to knowledge of Z. marina associated fungi and highlights a need for more investigation into the functional and evolutionary roles of microbial eukaryotes associated with seagrasses.

Author Correction: Metagenome-assembled genomes provide new insight into the microbial diversity of two thermal pools in Kamchatka, Russia.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Characterization of the Mycobiome of the Seagrass, Zostera marina, Reveals Putative Associations With Marine Chytrids.

Seagrasses are globally distributed marine flowering plants that are foundation species in coastal ecosystems. Seagrass beds play essential roles as habitats and hatcheries, in nutrient cycling, and in protecting the coastline from erosion. Although many studies have focused on seagrass ecology, only a limited number have investigated their associated fungi. In terrestrial systems, fungi can have beneficial and detrimental effects on plant fitness. However, not much is known about marine fungi and even less is known about seagrass associated fungi. Here we used culture-independent sequencing of the ribosomal internal transcribed spacer (ITS) region to characterize the taxonomic diversity of fungi associated with the seagrass, Zostera marina. We sampled from two Z. marina beds in Bodega Bay over three time points to investigate fungal diversity within and between plants. Our results indicate that there are many fungal taxa for which a taxonomic assignment cannot be made living on and inside Z. marina leaves, roots and rhizomes and that these plant tissues harbor distinct fungal communities. We also identified differences in the abundances of the orders, Glomerellales, Agaricales and Malasseziales, between seagrass tissues. The most prevalent ITS amplicon sequence variants (ASVs) associated with Z. marina tissues could not initially be confidently assigned to a fungal phylum, but shared significant sequence similarity with Chytridiomycota and Aphelidomycota. To obtain a more definitive taxonomic classification of the most abundant ASV associated with Z. marina leaves, we used PCR with one primer targeting a unique region of this ASV's ITS2 and a second primer targeting fungal 28S rRNA genes to amplify part of the 28S rRNA gene region corresponding to this ASV. Sequencing and phylogenetic analysis of the resulting partial 28S rRNA gene revealed that the organism that this ASV comes from is a member of Novel Clade SW-I in the order Lobulomycetales in the phylum Chytridiomycota. This clade includes known parasites of freshwater diatoms and algae and it is possible this chytrid is directly infecting Z. marina leaf tissues. This work highlights a need for further studies focusing on marine fungi and the potential importance of these understudied communities to the larger seagrass ecosystem.

Fungi in the Marine Environment: Open Questions and Unsolved Problems.

Terrestrial fungi play critical roles in nutrient cycling and food webs and can shape macroorganism communities as parasites and mutualists. Although estimates for the number of fungal species on the planet range from 1.5 to over 5 million, likely fewer than 10% of fungi have been identified so far. To date, a relatively small percentage of described species are associated with marine environments, with ∼1,100 species retrieved exclusively from the marine environment. Nevertheless, fungi have been found in nearly every marine habitat explored, from the surface of the ocean to kilometers below ocean sediments. Fungi are hypothesized to contribute to phytoplankton population cycles and the biological carbon pump and are active in the chemistry of marine sediments. Many fungi have been identified as commensals or pathogens of marine animals (e.g., corals and sponges), plants, and algae. Despite their varied roles, remarkably little is known about the diversity of this major branch of eukaryotic life in marine ecosystems or their ecological functions. This perspective emerges from a Marine Fungi Workshop held in May 2018 at the Marine Biological Laboratory in Woods Hole, MA. We present the state of knowledge as well as the multitude of open questions regarding the diversity and function of fungi in the marine biosphere and geochemical cycles.